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ABSTRACT 



43 

Detailed studies of resolved young massive star clusters are necessary to determine their dynamical state and evaluate the importance 
of gas expulsion and early cluster evolution. In an effort to gain insight into the dynamical state of the young massive cluster R 136 and 
obtain the first measurement of its velocity dispersion, we analyse multi-epoch spectroscopic data of the inner regions of 30 Doradus 
+jf ■ in the Large Magellanic Cloud (LMC) obtained as part of the VLT-FLAMES Tarantula Survey. Following a quantitative assessment 

of the variability, we use the radial velocities of non- variable sources to place an upper limit of 6kms~' on the line-of-sight velocity 
dispersion of stars within a projected distance of 5 pc from the centre of the cluster. After accounting for the contributions of undetected 
binaries and measurement errors through Monte Carlo simulations, we conclude that the true velocity dispersion is likely between 
4 and 5kms~' given a range of standard assumptions about the binary distribution. This result is consistent with what is expected 
if the cluster is in virial equilibrium, suggesting that gas expulsion has not altered its dynamics. We find that the velocity dispersion 
, would be ~25 km s~ ] if binaries were not identified and rejected, confirming the importance of the multi-epoch strategy and the risk 

£Nj ■ of interpreting velocity dispersion measurements of unresolved extragalactic young massive clusters. 
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1 . Introduction significant role in the early evolution/disruption of star clusters 

still needs to be addressed. 

The expre ssion "infant mortal ity" of star clusters was initially 

coined by [Lada & Lada| fl2003|) to describe the discrepancy be- Star clust ers have been observed to expand in fheir 



tween the number of observed open clusters and the number of first lOOMyr (iMackev & GilmoreL l200l iBastian et all [2008; 



1 embedded clusters. These authors argued that there are about ten |p or tegies Zwart et al 11201 Ol) . but this expansion is not direct ev 



times fewer open clusters than expected if all embedded clus- idence for the importance of gas expulsion. There are two ways 

ters evolve into open clusters. The rapid removal of gas left-over f or cl usters to expand as a response to mass loss (e.g. Hills, 

from star formation was suggested to explai n the apparent dis- |T980): (i) expansion following impulsive mass loss, e.g. change 

ruptioti of such a large fraction of clusters (e. g.lGever & BurkerA f potential due to removal of mass faster than the crossing time 

|200lt |Kroupa & Bojlyj |2002t |Bastian & Goodwirj |200fj|). The of the cluster, leaving the cluster in a super-virial state for a 

importance of the infant mortality scenario ho wever depends few crossing times, or (ii) adiabatic expansion, e.g. driven by 

on the definitio n adopted for embedded clusters (jBressert et alj, stellar evolution on a slow timescale compared to the crossing 

|2010t iBastianj |201l|), with more conservative criteria requir- time of the cluster (~10My r for young massive clusters, e.g. 

ing less than 50% of clusters to be destroyed to match the ob- |p or tegies Zwart et all l20T(l . in which case the cluster remains 

served number of open clusters. But no matter which definition in virial equilibrium. Thus, the best way to evaluate the impor- 

is adopted, the question of whether or not gas expulsion plays a tance of rapid gas expulsion (case i) and the implications for the 

formation and early evolution of star clusters is to determine the 

* Based on observations collected at the European Southern dynamical state of young clusters. In particular, it is important 

Observatory under program ID 182.D-0222 to verify if clusters are in virial equilibrium from a young age. 
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Attempts to determine the dynamical state of young clus- 
ters have been made by comparing dynamical masses (obtained 
through measuring the velocity dispersion and size of a cluster) 
and photometric masses (estimated from the age and integrated 
luminosity). For several unresolved extragalactic star clusters 
with ages of less than ~10Myr, the dynamical mass has been 
found to be up to ten times larger than the photometric mass (e.g. 
Basti an et al.L 12006). This led to the suggestion that these clus- 
ters might be super- virial and e xpanding following gas expul- 
sion dGoodwin & BastianLl2006l) . However. lGieles et al] d2010b) 
showed that the increase in the measured line-of-sight velocity 
dispersion in these young clusters could be produced by the or- 
bital motions of massive binaries. Massive stars indeed dominate 
the light of young massive clusters and their binary fraction is 
high (e.g. iMasonetali [2009; Bar ba et al.L l20Tot ISana & Evansl 
2011, Sana et al., submitted - hereafter Paper VIII). Therefore, 
the role of gas expulsion cannot be investigated through obser- 
vations of unresolved extragalactic clusters. 

Detailed dynamical studies of individual resolved clusters 
(i.e. in the Galaxy, Small Magellanic Cloud - SMC - or Large 
Magellanic Cloud - LMC) are necessary to make progress. In the 
case of radial velocity (RV) surveys, multi-epoch observations 
are needed to detect binaries, which can then be removed and 
allow a cleaner estimate of the velocity dispersion (unless the 
orbital solution is known, in which case the centre-of-mass ve- 
locity can be included in the dispersio n calculation). The y oung 
massive cluster R136 (M - 10 5 M n . lAndersen et all l2009t) . in 
the 30 Doradus region of the LMC, is an ideal target to test the 
impact of gas expulsion given its young age of less than 2Myr 
dde Koter et all Il998t iMassev & Hunter! Il998t ICrowther et all 
I2010I) . 

The identification of binaries was central to the study of 
stellar dynamics in 30 Doradus wi th the Gemini Multi-Object 
Spectrograph bv lBosch et al.l d2009l) . who found that the veloc- 
ity dispersion of NGC 2070 went from ~30km s _1 to 8.3 km s _1 
after correcting for orbital motions, confirming their prediction 
dBosch et all 1200 ll) that the high velocity dispersion of NGC 
2070 could be due to undetected binaries. In order to perform 
a similar study for the dense surroundings of R136, a different 
observational approach was required because this region is too 
crowded for fibre spectroscopy. A ke y componen t of the VLT- 
FLAMES Tarantula Survey (VFTS; lEvans et all 1201 ll here- 
after Paper I) is multi-epoch spectroscopy in the inner part of 
30 Dor with the FLAMES-ARGUS integral-field unit (IFU), 
which we use here to obtain an estimate of the velocity dis- 
persion of 'single' stars in R136. The central velocity disper- 
sion is an important proxy for the dynamical state of the clus- 
ter and is required to estimate, for example, the central poten- 
tial and relaxation time-scale, both of which are required for 
detai led numerical (TV-body) calculations of R136-li ke objects 
(e.g. iPortegies Zwart et~al l ll999falQieles et al.Ll2010ab . As a re- 
sult of the same study, the identification of massive binaries 
could provide important clues to the star-formation process and 
the subsequent dynamical evolution of R136, in which binaries 
can be formed, ejected and destroyed. Our measurements will 
also serve as useful empirical input for the modeling of stel- 
lar interactions in dense cluste rs, s uch as the recent stu dies of 
IFuiii & Portegies Zwartl d2011l) and lBaneriee et all d2012l) . 

We report here on a velocity dispersion estimate for the 
young massive cluster R136. We describe the FLAMES- 
ARGUS IFU data in Sect. [2] In Sect. [3] we present our RV 
and variability analysis of the stars observed with ARGUS. In 
Sect. [U we discuss the spectral classification of the non-variable 
ARGUS sources. In Sect. [5] we briefly introduce the VFTS 



Medusa data complementing the ARGUS data. We calculate the 
velocity dispersion from the stars showing no variability and also 
estimate the contribution of errors and undetected binaries in 
Sect. [6] The implications of the measured velocity dispersion 
for the dynamical state of R 136 are discussed in Sect. [7] Finally, 
we present our conclusions in Sect. [8] 

2. ARGUS data 

The main data used in this paper consist of five ARGUS IFU 
pointings in the central arcminute of 30 Dor (see Figure 2 of 
Paper I), for which at least five epochs were obtained. This cen- 
tral region is relatively gas free, so the nebular contamination in 
the spectra is not as important as in other regions of 30 Dor. 
Spectra of 41 different sources, corrected to the heliocentric 
frame, were extracted from the data cubes. To probe the stellar 
dynamics of R136, the LR02 setting of the Giraffe spectrograph 
was used (,1=3960-4570 A, AA = 0.40 A, R= A/AA -10500). 
This setting gives access to several stellar absorption lines suit- 
able for RV analysis (see section[3]) at the typical signal-to-noise 
ratio (S/N) of ~85 of our single-epoch spectra. To optimize the 
detection of binaries, the first two epochs were observed without 
time restrictions, and the third and fourth were observed with a 
minimum interval of 28 days between both the second and third, 
and third and fourth epochs. The fifth epoch was observed at 
least one year after the first epoch. Some observations were re- 
peated due to changes in conditions and other operational issues, 
providing additional epochs for three of the five pointings. Note 
that expanding the time baseline of the observations enhances 
the chances of detecting long-period binaries, but these binaries 
are the ones that have a smaller impact on the velocity dispersion 
enhancement. Further details on the ARGUS data, the reduction 
and the extraction procedure can be found in Paper I. 

For the analysis presented in this paper, individual exposures 
were considered to belong to the same epoch if their start time 
was separated by less than one hour. The multiple exposures 
composing a single epoch were then averaged using the errors 
(propagated throughout the reduction process) as weights and 
performing a 5 cr clip around the median to remove remaining 
cosmic features. The spectra from individual exposures had al- 
ready been normalised as part of the extraction procedure. The 
resulting epochs, their modified Julian date, and the correspond- 
ing ARGUS pointings and exposures are listed in Table lA.il 

An effort was made to extract spectra preferably for sources 
that appeared single by comparing the ARGUS data cubes with 
an ar chival HST Wide-Field Camera Three (WFC3) F555 W im- 
age dDe M archi et all 12.01 ll) and identifying matching sources 
(see Paper I). Out of 41 ARGUS sources, 23 are dominated by 
one bright source from the WFC3 image. For these 23 ARGUS 
sources, only stars at least a factor of ten fainter (in data counts) 
are visible in the WFC3 image within the region covered by the 
ARGUS spatial pixels from which the spectrum was extracted. 
Even if those fainter stars were contaminating the ARGUS spec- 
trum, it is doubtful that they would contribute significantly to 
the helium absorption lines used for our RV analysis (section 
13.21 ). The ARGUS spatial elements for a further 11 extracted 
sources are dominated by one bright object from the WFC3 im- 
age, but could suffer from more significant contamination from 
nearby stars (typically at a level of ~ 20%). The remaining seven 
ARGUS sources appear multiple in the WFC3 F555W image, 
with two or more densely-packed bright stars contributing at a 
comparable level to the flux in the region of the ARGUS source. 
These seven ARGUS sources were retained because they could 
still prove useful to our analysis (see below). 
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Note that even the spectrum of apparently single ARGUS 
sources (based on the WFC3 image) could contain contributions 
from multiple stars. The inner part of R136 is densely populated 
with stars. Two stars could be several thousand AU apart and, 
given the distance to the LMC, still appear as a single source 
in the WFC3 image. We would not expect to detect the mo- 
tion of unresolved binaries with such large separations given our 
spectral resolution and time coverage. These considerations are 
however not real concerns for our study because we can iden- 
tify shorter period binaries from RV variations (section l3~2l and 
estimate the residual contribution of undetected binaries to the 
velocity dispersion using Monte Carlo simulations (section lo3l l. 
Stars showing double/asymmetric line profiles or inconsistent 
absolute RVs between different lines can also be flagged as mul- 
tiple (either true binaries or a chance alignment of stars along 
the line of sight). On the other hand, we can still use the sources 
which appear multiple in the WFC3 image if they show none 
of the above spectroscopic signs of binarity/multiplicity. In that 
case, one star could be dominating the spectrum, or several stars 
could be contributing without showing any apparent RV differ- 
ence, and the source would still be valid to study the dynamics 
of the cluster. 



3. Radial velocity and variability analysis 

3.1. Zero-point errors 

To check that zero-point errors do not affect our multi-epoch 
RV measurements significantly, we first cross-correlated the 
ARGUS calibration arc spectra of each epoch/exposure with the 
arc spectrum of the corresponding fibre from the first exposure of 
the first epoch. We then determined the zero-point velocity shifts 
from the peak of the resulting cross-correlation functions. The 
corresponding distribution of zero-point errors for over 20 000 
such measurements is shown in Fig.Q](top panel). The vast ma- 
jority of the measurements lie between -0.5 and 0.5 km s with 
a peak at ~0 kms -1 , suggesting that the instrument wavelength 
calibration is remarkably stable and does not introduce spurious 
variations of stellar radial velocity between the different epochs. 
To investigate possible issues with the telescope and/or instru- 
ment not accounted for in the arc spectra, we also performed a 
similar analysis on the nebular lines (Hy and H<5) from selected 
spaxels with minimum stellar contamination. The resulting dis- 
tribution of velocity shifts for over 500 measurements peaks near 
km s _1 and has a small standard deviation of ~1 km s _1 (Fig.Q] 
bottom panel), which again suggests that there are no significant 
systematic shifts between epochs. 

3.2. Variability criteria 

In order to exclude the stars exhibiting RV variations from our 
calculation of the velocity dispersion, we implemented the fol- 
lowing quantitative assessment of the variability. 

For each star and epoch, we started by fitting Gaussian 
profiles to individual stellar absorption lines (Hen /i4200, 
Hen ,14542, Hei+Hen /i4026, Hei A4143, Hei ,14388 and 
He i ,14471, or a subset when some of these lines were too weak 
or the S/N too low). In principle, the profiles of the He lines are 
not Gaussian, particularly in the wings. However, given the mod- 
erate S/N of our spectra and the noise in the line wings, the use 
of Gaussians provided good fits to the line profiles (see Fig. [2} 
and did not affect the results. We used the MPFIT idl least- 
squares fitting routine ( MarkwardA l2009f) and Gaussian profiles 
defined by their central radial velocity, width, depth, and contin- 
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Fig. 1. Distribution of velocity shifts from the cross-correlation 
of ARGUS calibration arc spectra (top) and from the 
cross-correlation of nebular lines (bottom) between different 
epochs/exposures. 



uum value. An example fit is shown in Fig. [2] for the ARGUS 
source VFTS 1026. The resulting Icr error on the measured RV 
in this case is ±2.1 kms~', illustrating the precision that can be 
achieved for a single epoch with a good quality spectrum. To 
check if errors in the normalisation of the continuum could in- 
fluence the RV measurements, we also performed the fits allow- 
ing for a linear component instead of a constant continuum, but 
this did not affect the results. We chose the Gaussian fitting ap- 
proach (as opposed to cross-correlation for example) because it 
can provide reliable error estimates on the velocities and directly 
yields absolute RVs, which we need to compute the velocity dis- 
persion if different lines are used for different stars (see section 
13.31 ). This method is also well suited to our data set, for which 
the quality of spectra varies significantly between stars and even 
from one epoch to the next. 

The lines listed above are characteristic of the O-type stars 
composing the vast majority of our sample (see section [4] and 
Table [TJ. They were the only lines generally strong enough to 
perform satisfying RV measurements on the spectra of individ- 
ual epochs. Metallic lines were not considered because they were 
generally too weak, and Balmer lines were ignored because their 
profiles and apparent RVs can be affected by stellar winds or 
strong/variable nebular contamination. When one of the fitted 
lines appeared significantly contaminated by nebular emission it 
was also rejected, unless the nebular component could be clearly 
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Fig. 2. Example of a simultaneous Gaussian fit (i.e. same RV 
for all lines) to the Hen ,14200, Hei ,14388, and Hen ,14542 
absorption lines for an individual epoch of the ARGUS source 
VFTS 1026 (Q2-4.5 + mid/late O, see section!). 



identified and the fit could be performed on the wings of the line 
only. The most affected line was Hei ,14471, but it was also oc- 
casionally a problem for the other He i lines. 

From close inspection of all the fitted profiles, we first iden- 
tified the stars showing double-lined spectroscopic binary (SB2) 
profiles. Four of the 41 ARGUS sources show SB2 profiles 
(VFTS 1016, 1019, 1031 and 1033), and none of them had a 
primary component with a constant RV throughout all epochs 
(based on the criteria below), suggesting that they are genuine 
multiple systems and not just the result of the alignment of two 
stars with different RVs along the line of sight. The epochs at 
which SB2 profiles are observed in these sources are indicated 
in Table lA~2l 

We considered a star as a radial velocity variable if the se- 
ries of RV measurements of any of the fitted lines of that star 
contained at least one strongly deviant point, which we defined 

as 



RV~/i 



CTi 



>4, 



(1) 



where RVi and cr, are the radial velocity and its 1 cr error at epoch 
i, and // is the weighted mean RV over all the epochs. If the above 
condition was fulfilled for at least one individual line, then the 
null hypothesis of constant RV was rejected. 

For each star and each line, we also computed the value of 
X 2 assuming a constant RV, i.e. 



X 



-Z 



(RVi -/if 



(2) 



and rejected the constant RV hypothesis when the goodness of fit 
of a constant RV model to the data was poor, which we defined 

as 



1 -P(x 2 ,v) < 10~ 4 , 



(3) 



where P(x 2 , v) is the probability that, in a x 1 distribution with 
v degrees of freedom (v =# of epochs - 1), the value of x 2 is 
less than or equal to the value computed in Eq. [2] The thresh- 
olds adopted in Eqs. [TJ and [3] were chosen so that the proba- 
bility of a false variability detection in our sample (given the 
sampling and accuracy of our measurements) remained negligi- 
ble. In their analysis of the multiplicity properties of the O-type 



stars in VFTS, Sana et al. (Paper VIII) adopted slightly different 
variability criteria, but mention that their results are generally 
equivalent to those of a variability test based on the goodness of 
fit of a constant RV model like that of Eqs. [2] and [3] 

As a further check, we investigated line profile vari- 
ations (LPV) by computing Time Variance Spectra (TVS; 
I Fullerton et all 119961) . The only difference with the method of 
iFullerton et al.l ( 1 19961) is that we used the known error bars at 
each pixel to define a wavelength-dependent threshold instead 
of using a flat threshold based on continuum noise. This has the 
advantage of taking into account the varying S/N as a function 
of wavelength. The TVS at wavelength A is given by 



TVS (A) = 



N 



N-lZil/cr 2 (A) 



Z 



C/Kj) - if (A))) 2 



(4) 



where N is the number of spectra in the time series, fi(A) and 
o-j(A) are the flux and 1 cr error at wavelength A for epoch i, and 
(f(A)) is the weighted average flux of all epochs at wavelength 
A. The threshold for variability at a confidence level of 99% is 



TVS ( A) > 



N 



P y z(QM,N-l) 



Zil/cr 2 (A) 



N-l 



(5) 



where /^(O.Ol , N-l) is the cutoff value in &x distribution with 
N-l degrees of freedom such that the probability that a random 
variable is greater than this cutoff value is equal to 0.01. The 
cases for which significant variability is inferred from the TVS 
are indicated in Table lA.2l These results generally confirm those 
obtained from the other variability tests, and in two cases help 
to establish significant variability in emission-line stars with no 
or only weak absorption lines. Fig. [3] shows examples of TVS 
for two stars of our ARGUS sample, one revealing no variabil- 
ity and the other showing significant LPV. Note that the TVS 
often shows variability in nebular emission lines due to changes 
in conditions affecting the sky subtraction, but we ignore these 
spectral regions when assessing the variability of a star. The 
quantity TVS 1 ^ 2 is plotted in Fig. [3] because it scales linearly 
with the size of the spectral flux deviations and gives a direct 
estimate of the amplitude of the variations. 

If none of the three tests above (Eq. [TJ Eq. [3] TVS) re- 
vealed variability in individual lines, we performed simultane- 
ous Gaussian fits of He n ,14200, He n ,14542 and He i ,14388 (or 
a subset of these lines when one of them was too weak or sim- 
ply not present) by forcing their central RV to be the same (see 
Fig.[2]for an example). These three lines were adopted because 
they give reliable absolute RVs (see section [3~3l . We can there- 
fore fit them together to obtain more precise RV measurements. 
The series of RV measurements resulting from these simultane- 
ous fits were then tested for variability using Eqs.Q]and[3]again 
(see results in Table lA.2l . Stars still showing no sign of variabil- 
ity were then considered as suitable to study the dynamics of the 
cluster. The RVs of the ARGUS sources for all individual epochs 
are presented in Table lA.3l 

Out of 41 ARGUS sources, 16 are not detected as variable, 
17 are variable (four SB2, 11 SB1, two emission-line stars with 
variability determined from the TVS), seven have a too low 
S/N for a meaningful variability analysis, and another one is an 
emission-line star with no suitable absorption line for RV anal- 
ysis and for which no significant variability was detected from 
the TVS. 
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Fig. 3. Weighted mean spectrum (first and third panels from 
the top) and Temporal Variance Spectrum (TVS 1 '' 2 ; second and 
fourth panels from the top) for the ARGUS sources VFTS 1028 
and VFTS 1025. The red dashed curves indicate the 99% con- 
fidence level for variability. Significant variability is only seen 
in the nebular emission lines of VFTS 1028, but it is detected in 
several stellar lines in VFTS 1025. 



3.3. Absolute radial velocities 

We want to use as many lines as possible to increase the preci- 
sion of RV measurements, but at the same time we need to make 
sure that the selected lines give accurate results. In particular, all 
the lines fitted simultaneously to measure the RV of a given star 
should give consistent results when fitted individually, and the 
selected lines should provide consistent absolute RVs if differ- 
ent subsets of lines are used for different stars (e.g. He i lines for 
late O-type stars and He n lines for early O-type stars). 

Our choice of suitable lines for absolute RV measurements 
(He i ,14388, Hen ,14200, and Hen ,14542) is supported by the 
RV analysis performed in Paper VIII on the large sample of O- 
type stars observed with Medusa (see section |5j. The final RVs 
adopted for the non-variable ARGUS stars are obtained from 
simultaneous fits of these three lines or a subset of them. 

In stars with strong winds, wind infilling of photospheric 
lines can modify the RV measured by Gaussian fitting for lines 
like He n ,14200 and He n ,14542 and also result in RV shifts be- 
tween these lines. Comparison of Gaussian fitting measurements 
and values determined from CMFGEN models shows good con- 
sistency for O-type dwarfs and giants, but small shifts of a few 
tens of km s _1 for supergiants (P. Crowther, private communica- 
tion). In our calculation of the velocity dispersion (section loTTT i. 
we therefore payed particular attention to the stars that are iden- 
tified as possible supergiants (section|4]l. 



4. Spectral classification of ARGUS non-variable 
sources 

To get a general idea of the spectral content of our ARGUS sam- 
ple, we classified the non-variable and presumably single stars. 
We did not attempt the complex task of classifying the bina- 



ries/variable stars, partly because the limited wavelength cov- 
erage of the ARGUS spectra makes it even more difficult. 

To assign spectral types (SPT), we visually inspected the 
ARGUS spectra degraded to an eff ective resolving power of 
4000 and, following the premises of lSota et al.l <2011). we per- 
formed a morphological classification. In particular, for stars at 
intermediate and late subtypes, we used the eye-estimated line 
ratios of He i+n ,14026 to He n ,14200, He i ,1447 1 to He n ,14542, 
He i ,14388 to Hen ,14542, Hei ,14143 to Hen ,14200 and 
Si in ,14552 to He n ,14542 in order to assign the spectral subtype. 
For the hottest stars, on the other hand, we were not able to ex- 
ploit the primary criteria based on the Nitrogen ionization equi- 
librium because our ARGUS spectra do not cover wavelengths 
beyond ~4570 A. We thus concentrated on criteria related to 
the initial appearance of certain Hei lines (such as Hei ,14471, 
He i ,14143, and Hei ,14388), and the presence and strength of 
N rv ,14058 and the Si iv doublet in emission. The morphological 
classification of our targets was furthermore constrained by com- 
paring the degraded spectra to the spectra of O-type standards of 
solar metallicity compiled for the Tarantula Survey (Sana et al., 
in preparation) as well as to the spectra of VFTS targets obtained 
with Medusa-Giraffe which have already been classified. 

Assigning a luminosity class to the ARGUS targets was more 
difficult because we could not exploit the selective emission ef- 
fects in He n ,14686 and N m /l/i4634^1640-4642 (at subtypes 
earlier than 08), and the value of the He n ,14686/Hei ,14713 ra- 
tio (at subtypes 09-9.7). While at late-O types one could still 
rely on secondary criteria, such as the ratio of Si iv ,14089 to 
He i ,14026, no alternative luminosity diagnostics exist at early- 
and mid-O types. Given this situation, we decided to follow 
IConti & Alschulerl (119711) and exploit the increasing intensity of 
Si iv ,14089 relative to the nearby Hei ,14143 (at subtypes later 
than 05 only). 

The classification of the non-variable ARGUS targets, de- 
rived as outlined above, is presented in Table \T\ In section [B] 
of the appendix, we comment on specific sources, in particular 
those that have revised spectral types and the three that appear 
to have composite spectra. The accuracy of the spectral types 
reported in Table Q]for the ARGUS sources is typically between 
one and one and a half subtypes, with u ncertainties caused b y the 
effects of n ebular emission , rotati on dMarkova et al.L 1201 lb . and 
metallicity dMarkova et al.L 120091) . The uncertainty on the lumi- 
nosity class is significantly larger as we were only able to sepa- 
rate the stars into two broad categories: high luminosity objects 
(luminosity class I/II) and low luminosity objects (luminosity 
class III/V). For completeness, we also include previously pub- 
lished spectral types in Table [T] 



5. Supplementary Medusa data 

To complement our RV measurements of the ARGUS sources, 
we also include the RV measurements of the non-variable O- 
type stars observed with Medusa-Giraffe (Paper VIII) in the in- 
ner 10 pc of R136 . This gives us 22 additional stars, 20 of which 
are located between 5 and lOpc from the centre of the cluster. 
We do not consider stars beyond lOpc from the centre in our 
analysis of R136. Although somewhat arbitrary, this cutoff at 
10 pc is a reasonable trade-off between increasing the number of 
stars in our sample and limiting the possible contamination from 
nearby clusters or other star formation events in the surroundings 
ofR136. 

The Medusa observations (see Paper I) were performed 
using three of the standard Medusa-Giraffe settings (LR02: 
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Fig. 4. Distribution of ARGUS and Medusa sources used in this work overlaid on an F555 W HST-WFC3 image. ARGUS stars in 
which no variability was detected are shown with blue circles, ARGUS variable stars are represented by red squares, and ARGUS 
stars with too low S/N for RV analysis or no suitable absorption lines are indicated by red crosses. The ARGUS IFU pointings 
(Al to A5) are also shown as grey transparent rectangles. Medusa non-variable stars added to the sample are represented by green 
diamonds. The dashed circles indicate projected radial distances of 2, 4, 6, and 8 pc from the centre of R136. 



i=3960^564 A, AA = 0.61 A, R-7 000; LR03: ,1=4499- 
5071 A, AA = 0.56 A, R-8 500; HR15N: /t=6442-6817 A, 
AA = 0.41 A, R-16000). To detect RV variables, six epochs 
were observed for the LR02 setting with time constraints similar 
to those of the ARGUS observations (section|2]and Table lA.ll . 

The RV and variability analysis of these stars is presented in 
Paper VIII in the context of a study of the multiplicity of O-type 
stars in 30 Dor. The method is similar to the one we applied in 
section [3~2l As a consistency check, we applied our method on 
the LR02 Medusa spectra of several O-type stars and found ab- 
solute RVs and 1 <x uncertainties fully consistent with the results 
from Paper VIII. Thus, we do not repeat all the RV measure- 
ments of the Medusa stars, but adopt the values of Paper VIII 
instead. For two of those additional Medusa stars, the luminos- 
ity class could not be constrained, but none of the others is clas- 
sified as a supergiant (Walborn et al., in preparation). The RVs 
of these stars, obtained from Gaussian fitting, should therefore 
be reliable as discussed in section 1331 The spectral types of the 



non-variable Medusa O-type stars, determined by Walborn et al. 
(in preparation), are listed in TableQ] 

Note that in parallel to the ARGUS observations, a few 
stars were also observed in the inner lOpc of R136 with the 
Ultraviolet and Visual Echelle Spectrograph (UVES), provid- 
ing a greater resolving power than the Giraffe spectrograph. All 
the UVES targets but one in this region were also observed 
with ARGUS or Medusa, and these were already shown to 
be binaries/variable based on the ARGUS or Medusa observa- 
tions alone. One star was observed only w i th UV ES; Mk39, an 
02.5 If*/WN6 star dCrowther & WalbornL l201lh . which is not 
suitable for our RV analysis because the He n lines are wind con- 
taminated. 

Our final sample of apparently single ARGUS and Medusa 
stars within lOpc from the centre of R136 is presented in 
Table [T] These are all the stars that we use for the analysis of 
the dynamics in the following sections. Their absolute RVs and 
their projected distance from the centre of the cluster (which 
we adopt to be the position R136-al) are also listed. The pro- 
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Table 1. RVs of ARGUS and MEDUSA sources showing no significant variability. 



ID 




Data 


a 


6 


r d 


RV 


Spectral type 


VFTS 


Aliases* 




(2000) 


(2000) 


[pc] 


n —It 

[km s '] 




1024 


S99-147, R136-226, MH-623 


ARGUS 


05 38 42.685 


-69 06 07.03 


1.1 


263.7±2.3 


08.5 III/V 


1014 


S99-56, R136-29, MH-203, P93-863 


ARGUS 


05 38 41.515 


-69 06 00.83 


1.2 


266.5+0.7 


02-4.5+mid/late O 
03 V [MH98] 


1012 


S99-249, R136-151, MH-178 


ARGUS 


05 38 41.386 


-69 06 02.49 


1.3 


264.6+2.8 


09 III/V 


1009 


S99-165, R136-88, MH-141 


ARGUS 


05 38 41.163 


-69 06 02.83 


1.6 


274.2±1.4 


06.5 III/V 


1007 


S99-95, R136-60, MH-129, P93-827 


ARGUS 


05 38 41.077 


-69 06 01.74 


1.7 


264.1 ±1.3 


06.5 III/V 


1004 


S99-193, R136-126, MH-95 


ARGUS 


05 38 40.848 


-69 06 04.51 


2.0 


274.4±2.0 


09.5 III/V 


1008 


S99-163, R136-96, MH-134 


ARGUS 


05 38 41.108 


-69 05 58.33 


2.0 


277.3±1.0 


ON6.5 I/II 


1023 


S99-142, R136-83, MH-591 


ARGUS 


05 38 42.631 


-69 06 10.91 


2.0 


266.4±2.0 


08 III/V 

06 V [MH98] 


1026 


(S99-191), R136-41, MH-716, Mk35N 


ARGUS 


05 38 43.083 


-69 06 11.26 


2.2 


265.0±0.9 


02-4.5+mid/late O 
03 III(f*) [MH98] 
08: V [WB97] 


1002 


S99-312, R136-194, MH-77 


ARGUS 


05 38 40.646 


-69 06 05.73 


2.4 


272.7±5.6 


09.5 III/V 


1018 


S99-88, R136-37, MH-290, P93-900 


ARGUS 


05 38 41.887 


-69 06 12.45 


2.4 


258.4±0.6 


02-4.5+mid/late O 
03 III(f*) [MH98] 


1020 


S99-178, R136-101, MH-314 


ARGUS 


05 38 42.023 


-69 06 16.75 


3.4 


267.4±0.9 


03-4 


1006 


S99-257 


ARGUS 


05 38 41.066 


-69 06 15.16 


3.4 


264.2±3.9 


06.5 III/V 


1028 


S99-37, R136-23, P93-1036, Mk35S 


ARGUS 


05 38 43.274 


-69 06 16.45 


3.5 


271.9±0.5 


03.5-4.5 

03 III(f*) [MH98] 
04-5 V: [WB97] 


1010 


S99-189, R136-103, MH-156 


ARGUS 


05 38 41.268 


-69 06 16.94 


3.7 


283.1+2.2 


07 III/V 


468 


S99-86, MH-17, P93-706, Mk36 


Medusa 


05 38 39.38 


-69 06 06.39 


4.0 


272.0+0.8 


02 V((f*)) + OB 


477 


S99-455 


Medusa 


05 38 39.75 


-69 05 50.55 


4.5 


265.3±3.5 


0((n)) 


1035 


S99-169, R136-109 


ARGUS 


05 38 44.321 


-69 05 45.05 


5.0 


268.5±2.2 


08.5 I/II 

08 V [MH98] 


601 


S99-91, MH-986, P93-1317, Mkl4N 


Medusa 


05 38 46.29 


-69 05 59.25 


5.1 


266.8±0.6 


05-6 V((n))z 
05 V((f)) [MH98] 
04 V [WB97] 


484 


S99-124 


Medusa 


05 38 40.37 


-69 05 43.72 


5.3 


284.2±1.0 


06-7 V((n)) 


611 


S99-270 


Medusa 


05 38 46.90 


-69 05 58.71 


5.9 


265.4±2.2 


08 V(n) 


554 


S99-343 


Medusa 


05 38 43.79 


-69 05 38.70 


6.1 


275.6±1.4 


09.7 V 


446 


S99-194 


Medusa 


05 38 38.26 


-69 06 17.29 


6.4 


258.9±5.0 


O Vnn((f)) 


607 


S99-294 


Medusa 


05 38 46.76 


-69 05 48.75 


6.6 


257.8±1.0 


09.7 III 


476 


S99-206 


Medusa 


05 38 39.75 


-69 05 39.21 


6.7 


270.1±2.2 


0((n)) 


505 


S99-265 


Medusa 


05 38 41.49 


-69 05 34.52 


7.0 


265.5±3.4 


09.5 V-III 


536 


S99-295 


Medusa 


05 38 42.82 


-69 06 32.43 


7.2 


248.4±1.4 


06 Vz 


582 


S99-414 


Medusa 


05 38 45.19 


-69 05 37.43 


7.2 


270.2±2.3 


09.5 V((n)) 


615 


S99-218 


Medusa 


05 38 47.33 


-69 06 17.70 


7.3 


256.4±7.6 


09.5 Illnn 


515 


S99-434 


Medusa 


05 38 41.86 


-69 05 32.40 


7.4 


267.4±2.0 


06-9p 


513 


S99-266 


Medusa 


05 38 41.81 


-69 05 31.91 


7.6 


266.3+1.3 


06-7 Ill-II(f) 


622 


S99-333 


Medusa 


05 38 48.12 


-69 06 12.27 


7.8 


272.0±1.6 


09.7 III 


498 


S99-347 


Medusa 


05 38 41.19 


-69 06 35.17 


8.0 


265.9+7.2 


09.5 V 


451 


S99-346 


Medusa 


05 38 38.52 


-69 06 29.92 


8.3 


277.3+ 10.8 


09: ni:(n) 


483 


S99-309 


Medusa 


05 38 40.24 


-69 05 30.78 


8.3 


275.4+2.7 


09 V 


540 


S99-372 


Medusa 


05 38 43.08 


-69 06 36.88 


8.3 


242.7±1.1 


B0 V 


560 


S99-350, P93-1139 


Medusa 


05 38 44.35 


-69 06 40.73 


9.5 


259.7±0.9 


09.5 V 


465 


S99-365, P93-700 


Medusa 


05 38 39.28 


-69 06 38.93 


9.6 


255.8+4.4 


05: Vn 



References. MH = 



(S99-), Hun ter et al 



Massev & Hunted dl998l). WB97 = [Walborn & Blades) dl997h . 'Aliase s identification numbers are from lSelman et alJll999l 
19971 (R136-). lMaiumufh & Heap 1994 (MH-), Parker 1993 (P93-) and Melnick 1985 (Mk). 



Notes. Spectral types were determined in this work (see section[4]l for the ARGUS sources and by Walborn et al. (in preparation) for the Medusa 
sources, unless otherwise indicated. The sources are sorted by increasing projected radial distance (r d ) from R136-al (a = 5 h 38 m 42 s . 39, 6 = 
-69°06'02."91, J2000). 



jected distances are for an adopted distance modulus of 18.5 mag 
(see Paper I). The absolute RV adopted for a given star is the 
weighted mean RV of all epochs. The spatial distribution of vari- 
able and non- variable ARGUS sources is shown in Fig. [4] along 
with the positions of the non-variable Medusa O-type stars added 
to our sample. 



6. Velocity dispersion 

6. 1 . Velocity dispersion upper limit 

We can now determine the observed line-of-sight velocity dis- 
persion (o"i D ) from the RV measurements of the non-variable 
stars. Because any effect not yet taken into account (e.g. unde- 
tected binaries, intrinsic errors) would tend to increase the in- 
ferred cr 1D , we can consider that the results shown in this subsec- 
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Fig. 5. Observed line-of-sight velocity dispersion, as a function of R, for the stars within a projected radial distance R from the centre 
of R136. In the left panels, stars of all subtypes are included, while in the right panels, we divide them into two subsamples: 07 
and earlier (blue triangles), 08 and later (red squares). In the top panels, all the non-variable stars (see TableQ} are included. In the 
middle panels, two stars with a radial velocity more than 3<x away from the mean RV are excluded. In the bottom panel, supergiant 
and SB2 candidates are also excluded. In each panel (and subsample), the first point from the left is the velocity dispersion of the 
four stars (of that subsample) closest to the centre, the second point is the velocity dispersion of the innermost five stars, and so on. 



tion represent an upper limit to the actual line-of-sight velocity 
dispersion of the cluster. 

In what follows , crip is determin ed by computing the stan- 
dard deviation, i.e. (Bevingto^,[l96§ 



H = 



1 y RVj 



Var=^Y^!, 
y,l/cr 2 ^ a 2 



o"id = -i/Var 



N 



N-V 



(6) 



(7) 



(8) 



where <r, is the uncertainty on the RV measurement RVj, fi is the 
weighted mean RV, Var is the variance, and N is the number of 
measurements. 

In the upper left panel of Fig. [5] we present, as a function 
of projected radial distance from the centre of R136, the <x 1D of 
all the non-variable stars within that radius. Apart from an ap- 
parent increase in <x 1D in the inner region (most likely due to 
the low number of stars and associated large uncertainty on <x 1D 
and its error), the velocity dispersion profile is relatively flat. For 
the stars within 5pc from the centre, we find cr 1D < 6 kms -1 . 
Two stars dominate the increase in cr 1D between 5 and lOpc: 



VFTS 536 (r d =7.2pc; RV=248.4±1.4 kms" 1 ) and VFTS 540 
(r d =8.3pc; RV=242.7+1.1 kms- 1 ). 

As a next step, to see the effect of possible outliers (slow run- 
aways or massive stars along the line of sight but not members 
of R136), we exclude the stars with an RV more than 3 cr away 
from the (weighted) mean RV of our sample (267.7 kms -1 ). We 
choose cr to be 6 kms -1 , the observed dispersion of the stars 
within 5pc from the centre. VFTS 536 and VFTS 540 are in- 
deed >3cr outliers, and excluding them results in an even flatter 
profile (Fig. [5] middle left panel). 

When also excluding possible supergiants (I/II or no lumi- 
nosity class attributed in Table Q]), for which the RVs obtained 
from Gaussian fits could be problematic (see section 13.3b . and 
SB2 candidates (composite spectra, see Table [TJ, the results do 
not change significantly, although more fluctuations are seen in 
the profile and the error bars are larger due to the smaller num- 
ber of stars (Fig. [5] bottom left panel). The apparent increase in 
the inner 2-3 pc also disappears when these supergiant and SB2 
candidates are excluded. 



In section 1631 we estimate the contribution of errors and un- 
detected binaries to cr 1D and attempt to reproduce the observed 
velocity dispersion for the stars within 5 pc (in projection) from 
the centre, i.e. cr 1D = 6 kms 1 . We could choose a different ra- 
dius, but because the velocity dispersion profile appears remark- 
ably flat in the inner lOpc, this would not change the results 
significantly. There is also a natural cut at ~5 pc if we consider 
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the de finition of a cluster proposed bv lGieles & Portegies Zwartl 
(1201 ll) . which states that stellar agglomerates for which the age 
of the stars exceeds the crossing time are bound and thus referred 
to as star clusters. The crossing time (i.e. the distance for a star 
to travel from one side of the cluster to the other; 2r) is roughly 
<t 1d x age. Given an age of ~2Myr and <x 1D ~ 6 kms -1 , we can 
conclude that the stars that are physically within ~6 pc from the 
centre are part of the cluster following the above definition. 



6.2. Contamination by "halo" stars 

The light profile of R136 suggests that it is not a single- 
component cluster but the composite of a real cluster and a 
"halo", i.e. an OB association, with the latter contributing to 
more than 5 0% and possibly as much as 90% of its total inte - 
grated light dMafz-A pellaniz, 1200U iMackev & Gilmorel 12003). 
It is therefore worth asking how much could the OB association 
contaminate our velocity dispersio n measurement fo r the cluster. 
In the double-comp onent EFF fit (lElson et aT , 1 19871) to the light 
profile of R 136 by Mac kev & Gilm ore (2003), the projected ra- 
dius where the two components contribute equally is at about 
5 pc. If we consider the inner profile as the cluster and the outer 
profile as the halo, we can conclude from this fit that the contri- 
bution of halo stars projected on the cluster is negligible (< 5%) 
in the inner 1 .25 pc, but it could be significant beyond 5 pc. 

We might expect the OB association to have a low velocity 
dispersion, as is observed in Galactic cases once binaries and 
runaways are excluded (e.g. lde Bruiinel[T999l) . However, the lat- 
ter are hard to identify if there is a massive cluster with a higher 
velocity dispersion embedded in the OB association. To explore 
the cluster/halo dichotomy, we divided our sample into two sub- 
type groups such that they had roughly the same number of stars, 
which resulted in these two subsamples: earlier than 07 and later 
than 08. We then computed the velocity dispersion as in sec- 
tion [6J] for the stars that could be placed in one of these groups 
(Fig. [5] right panels). We would expect the earlier type stars to 
be more concentrated towards the centre as a result of mass seg- 
regation or due to an overall age difference between the cluster 
and halo populations. Fig.|5]indeed suggest a higher concentra- 
tion of early-types towards the core, but this could well be due to 
the increasing crowding effects in the innermost regions favor- 
ing the detection of the brightest (most probably earliest) stars. 
In any case, there is no obvious difference in velocity dispersion 
between the two subsamples, and the contribution of halo stars 
to the velocity dispersion remains very difficult to identify. 



6.3. The contribution of cluster rotation 



In an accompanying letter dHenault-Brunet etal we 
present evidence for rotation of R136. From comparison of our 
RV measurements with different simple rotating models, we in- 
fer a rotational amplitude of ~3kms~ 1 and an optimal position 
angle for the rotation axis at an angle of ~45° east of north. To re- 
move the anisotropy due to the suggested rotation and its contri- 
bution to the computed velocity dispersion, we subtracted from 
our measured RVs the rotation curve from these simple models. 
We find that the velocity dispersion obtained after this correction 
is typically 0.5kms~' lower than the values presented in sec- 
tion [6T| Thus, a small component of the observed line-of-sight 
velocity dispersion could be attributed to cluster rotation. 



6.4. The velocity dispersion when including binaries/variable 
stars 

It is interesting to see what the computed ctid would have been if 
we had not been careful about identifying and rejecting binaries 
and variable stars, or if we were not dealing with multi-epoch 
observations. 

We randomly selected one epoch for each ARGUS source 
on which the RV analysis was performed (variable and non- 
variable sources), and repeated the process for 10000 combina- 
tions of epochs. A median <x ID of 25.0kms _1 was obtained from 
all these combinations of single-epoch RV measurements, with a 
standard deviation of 5.9 km s _I and values of cr 1D ranging from 
12.9 to 48.0km s -1 . If we do a similar test but limit ourselves 
to the non-variable ARGUS sources (Table[TJ, we find a median 
<x ID of 6.2kms~ 1 (in good agreement with the results of section 
16. U . with a standard deviation of 0.7 km s _1 . 

Note however that the velocity dispersion obtained when in- 
cluding all the RV variables cannot be directly compared with 
the velocity dispersion one would obtain from the integrated 
light of a distant star cluster, as a few outliers could increase 
the velocity dispersion significantly without contributing much 
to the integrated light. 

6.5. The contribution of errors and undetected binaries 

To estimate the contribution of measurement errors and un- 
detected binaries to the observed velocity dispersion, we per- 
formed a series of Monte Carl o simulations. Th ese are adapted 
from the method presented by ISana et all (|2009) and refined in 
Paper VIII to e s tima te the probability to detect binary systems. 
Bosch & Meza (l200ll) have also previously performed similar 
Monte Carlo simulations testing intrinsic and observed proper- 
ties of binary stars. Our simulations mimic the process that we 
have been going through, i.e. identify variables from series of RV 
measurements and then compute the velocity dispersion from the 
remaining non-variable stars. 

The general procedure goes as follows. We first adopt rea- 
sonable orbital parameter distributions (period, mass ratio, ec- 
centricity) and an intrinsic binary fraction. Then, for 10 000 pop- 
ulations of N stars (where N is the size of the sample, i.e. all the 
ARGUS sources plus Medusa O-type stars within 5 pc from the 
centre of R136 on which the RV variability analysis was per- 
formed), we randomly draw which are binaries and which are 
single, and also randomly draw the parameters for the binaries 
from the adopted distributions. From these, we compute the or- 
bital velocity at each epoch (based on the time sampling of the 
observations) for all the binaries assuming random orientations 
of the orbital planes and uncorrelated random time of perias- 
tron passage. We then add the measurement noise (based on the 
RV uncertainties of our sample) to the computed RVs of bina- 
ries and single stars. By applying the RV variability criteria of 
section [3~2l we can then eliminate the stars that we would have 
flagged as variable, compute the line-of-sight velocity dispersion 
of the apparently non-variable stars, and estimate the contribu- 
tion from the orbital motion of binaries that were not detected. 
This procedure has the advantage that we not only take into ac- 
count the long period binaries (i.e. too long a period to be de- 
tected with the VFTS), but also all the shorter binary systems 
missed by our time sampling, with a statistical weight exactly 
defined by the incompleteness. 

Before performing the full procedure outlined above, we es- 
timated the contribution of measurement errors by taking a pop- 
ulation with no binaries. For each star, we drew the RV at each 
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Fig. 6. Estimate of the line-of-sight velocity dispersion distribu- 
tion for massive stars in the inner 5pc of R136. The blue dot- 
ted curve shows the contribution of measurement errors. The red 
dashed curve is the distribution from binaries undetected after 
applying our variability criteria (the initial input population has 
a binary fraction of 100%, and an Opik-law distribution of the 
periods is adopted with period range of 0. 15-6.85 in log P with 
P in units of days). The green dashed-dotted curve is the dynam- 
ical velocity dispersion which best reproduces the observed ve- 
locity dispersion. The black solid curve takes into account mea- 
surement errors, undetected binaries, and the dynamical velocity 
dispersion. The median (central tick) and 68% confidence inter- 
val (equivalent to +1 cr for Gaussian distributions) of the distri- 
butions are indicated on the upper part of the graph. 



epoch from a Gaussian distribution centered on zero and a sigma 
corresponding to the RV error at the corresponding epoch. We 
then applied the RV variability criteria of section 13.21 to make 
sure that our adopted thresholds do not lead to false detections 
(the false detection rate was indeed found to be negligible). We 
finally computed the velocity dispersion of the population and 
repeated this for 10000 populations. The resulting velocity dis- 
persion distribution is represented by the blue dotted curve in 
Fig. [6] The peak of this blue curve is at 0.80 (kms -1 ) -1 , out of 
the graph. It shows that, given the precision of our RV measure- 
ments, the intrinsic contribution of errors to the observed veloc- 
ity dispersion is small. Note that we performed two tests on the 
population size: (i) we first used N stars, the full population size, 
and (ii) we randomly picked stars from the initial population fol- 
lowing the results of a binomial test with a success rate of 50%, 
mimicking the fact that the dispersion is usually computed using 
an effective population of about half the original one because of 
binary rejection. The two approaches made no difference on the 
resulting contribution of measurement errors. 

The contribution of undetected binaries depends on what is 
assumed for the input distributions. We focus here on the distri- 
bution of orbital periods, because the distributions of mass ratios 
and eccentricities have a limited impact (see Paper VIII). In what 
follows we adopt a fla t distribution of mass r atios, as motivate d 
by t he recent results oflSana & Evansl (1201 lh . lSana et all d 201 2|) . 
and Kiminki & Kobulnickvl d20 1 2lK 

As a first test to estimate the intrinsic contribution of un- 
detected binaries, we adopted a conservative binary fraction of 
100% and a standard Opik-law distribution for the period (i.e. 
flat distribution of logP) with a period range of 0.15-6.85 (in 
log P, where P is in units of days). The maximum period adopted 



corresponds to the extrapolation of the Opik law until a 100% 
binary fraction is reached when considering the observed binary 
fraction and the overall detection rate of VFTS (Paper VIII). We 
ran the full procedure outlined above, but without adding the 
measurement noise to the extracted orbital velocities. The re- 
sulting velocity dispersion distribution (from the stars that are 
not identified as RV variable) is shown with the red dashed curve 
in Fig. [6] Under these assumptions, the velocity dispersion from 
undetected binaries is 4.6+jj km s~' . The quoted value is the me- 
dian of 10000 populations, and uncertainties correspond to the 
68% confidence interval (equivalent to +1 cr for Gaussian distri- 
butions). Note that within uncertainties, the undetected binaries 
alone can produce the observed velocity dispersion of 6 km s . 

To recover the true velocity dispersion of the cluster, we re- 
peated the simulation, this time adding the measurement noise 
to the orbital RVs. We also included a contribution from the dy- 
namical velocity dispersion that was varied until the most prob- 
able velocity dispersion matched the observed velocity disper- 
sion. For simplicity, we assumed that this dynamical contribu- 
tion to the velocity dispersion did not vary as a function of ra- 
dius. We estimate that the line-of-sight velocity dispersion at- 
tributable to the cluster dynamics is 3.9*q5 kms _I . The errors 
on this value correspond to changes in the input dynamical ve- 
locity dispersion that result in median values of the output ve- 
locity dispersion distribution at percentiles 0.16 and 0.84 of the 
optimal output distribution (i.e. the 'overall' simulated distribu- 
tion for which the median is the observed velocity dispersion). 
Note that we obtain a velocity dispersion of 26+9 km s~' from 
these simulations for a single epoch when including all binaries 
(i.e. without applying our variability criteria and thus without 
rejecting the binaries that would be detected), in keeping with 
the value of ~25kms~ 1 that we found in section l6~4l from the 
single-epoch RV measurements of all ARGUS sources (variable 
and non-variable). 

We also ran simulations using different binary fractions and 
period distributions. We first considered only periods shorter 
than 10 3 5 days, i.e. the ones that could be detected by VFTS 
(Paper VIII). Assuming a 50% binary fraction and a standard 
Opik law for the period distribution, the estimated velocity dis- 
persion from the undetected shorter period binaries alone is 
2.4^[ q kms~'. If instead of the Opik law we adopt the distri- 
bution measured in Paper VIII for the VFTS O-type binaries 
(/(logP) oc (logP)" 045 ), we obtain a velocity dispersion from 
the undetected shorter period binaries of 2.0^ ^ kms _I . This is 
the minimal contribution of undetected binaries. These values 
(to be compared with the red dashed curve in Fig.|6]l indicate 
that the contribution of undetected binaries is dominated by long 
period systems outside the sensitivity range of the VFTS. 

Recall that the velocity dispersion from undetected binaries 
when including very long period systems as well (i.e. binary 
fraction 100%, Opik-law distribution of periods ranging from 
0.15 to 6.85 in logP) is 4.6+}'j , kms~ I . If we adopt the distri- 
bution with /(log P) oc (log p)~ - 45 instead, the velocity disper- 
sion from both short and long period undetected binary systems 
is 3.4^[ jkms"'. In that case, the dynamical velocity dispersion 
that best reproduces the observed velocity dispersion would be 
5 km s _1 . The Opik law might therefore overestimate the contri- 
bution of binaries by ~1 kms -1 compared to the period distri- 
bution measured in VFTS O-type binaries (Paper VIII), but we 
should bear in mind that the extrapolation of the period distribu- 
tion to long periods is very uncertain. 

In summary, we estimate that the velocity dispersion due 
to cluster dynamics alone (i.e. removing the effect of bina- 
rity but still including a small contribution from rotation) is 
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likely somewhere between 4 and 5 km s 1 , with a contribution 
of ~ 0.5kms~' from rotation (see section [63}. 



7. Discussion 

Now that we have our measurement of the velocity dispersion at 
hand we can test the hypothesis that the cluster is in virial equi- 
librium. This is often done by deriving a dynamical, or virial 
mass, from the velocity dispersion which is then compared to 
the photometrically determined mass. To be able to derive the 
former we require knowledge about the mass distribution of the 
stars. For equilibrium models with an isotropic velocity disper- 
sion the one-dimensional velocity dispersion, cr 1D , relates to the 
mass, M, and the viral radius, r v , as M = 6r v cr 2 D /G. The virial 
radius is defined as r v = GM 2 /(2W), with W the potential en- 
ergy of the system. This relation is often expressed in terms of 
the radius containing half the light in projection, or effective ra- 
dius (r e ff), as M — T]r e ffcr 2 D /G, with 77 w 10. This is under the 
assumption that light traces mass, that the half- mass radius (?y ) 
in projection is 3/4 times the 3D half-mass radius (Spitzer, 1987) 
and that the ratio r v /rh ~ 5/4. The first tw o assumptions are not 
valid when a cluster is mass segregated dFlecketall [2006) as 
the 2D light radius can be twice as sma ll as the 3D mass radius 
dGaburov & GielesLl2008HHurlevLl2007l) . For models with very 
fiat density profiles, it is difficult to estimate the half-light ra- 
dius. The surface brightness profiles of young clusters are often 
approximated by cored templates with a power-law decline of 
the form I(r) = 7(0)(1 + r 2 //^) ~ 7 ' 2 , where ro is a scale radius . 
These profiles are often referred to as EFF profiles dElson et al.L 
119871) ." For y > 2 these models contain a finite amount of light, 
but diverge to infinite luminosity when y < 2. The boundary 
at y = 2 corresponds to 3D light profiles that decline as r~ 3 . 
For y larger than, but close to 2, the ratio r v /rh becomes ver y 
sensitive to the exact value of y dPortegies Zwart et all 12010). 
Additionally, determining r e ff becomes difficult as this quantity, 
and the total luminosity, can become unrealistically large when 
extrapolating to infinity. 

The light profile of R136 has a profile close to y — 2 
dMcL aughlin & van der MareJ 120051) . with indi cations for a 
'bum p ' in the optical l i ght profile at about lOpc ( Hunter et all 
119951; iMafz-Apellanizl |200T1: iMackev & GilmoreL 120031) . The 
presence of an additional component from the larger scale and 
near-constant-density OB association in which R136 is located 
(see section 16.21 ) has been interpreted as the reason for this 
relatively flat profile. In the near-infra red (NIR), the profiles 
are even flatter than the critical value dAndersen et all l2009t 
Campbell et alll2010l) and the halo structure is not that obvious, 
but the NIR data did not extend very far into the OB association. 
Mass segregation, age differences between the core and halo, and 
differential extinction (which becomes important ~10pc away 
from the centre of R136) could explain the flatter profile in the 
NIR compared to the optical. 

With the above caveats in mind, it is still interesting to 
see what dynamical mass we obtain for R136. If we assume 
77 ss 10, 3.4 < cr 1D < 6.0 km s (section J6}, an d adopt a 
half-light radius of r e ff = 1.7 pc {Hunter et all 11995) which is 
consistent with the half-light radius obtained for the inner com- 
ponent of the double EFF fit discussed in section 16.21 we get 
M = 4.6 - 14.2 x 1O 4 M . This is c onsistent with the estim ated 
photometric mass of ~ 10 5 M a bv lAndersen et all d2009l) . for 
which the cluster mass of ~ 5 x 10 4 M (computed for stellar 
masses between 25 M and down to 2.1 M ) was extrapolated 
assuming a Salpeter slope down to 0.5 M . 



However, because of the difficulties outlined above, we de- 
cide to address whether R136 is in virial equilibrium by explor- 
ing an alternative method for which we do not need to know the 
total mass. This relies on estimating the central velocity disper- 
sion by following a very sim ilar method to that presented by 
iRichstone & Tremaind d!986l) . From the observations we find 
that the velocity dispersion is roughly constant with radius (sec- 
tion !6. 11 >. We can express the expected central dispersion ctid(0) 
of (self-consistent and isotropic) models with isothermal inner 
parts in terms of observed properties of the cluster 



cr 2 lD (0) = anGr v I(0)r , 



(9) 



where G « 0.0043 pcM km s is the gravitational constant, 
Yy the mass-to-light ratio in M /L in the V-band and a de- 
pends on the model. 

We first consider the modified Hubble profile dRood et all 
119721) . which is an EFF profile with y — 2, very close to the best 
fit to the surface brightness profile in the optical. From solving 
Jeans' equations assuming hydrostatic equilibrium and isotropic 
velocities we find that a = 3 - 41n(2) « 0.227. Secondly, we 
look at the lPlummerldl911l) model for which a = 1/6 * 0.167. 
Finally, we can consider the isothermal sphere, which cannot be 
expressed in an EFF profile. For this model ro is defined a£] 
r? = 9crj D /[47rGp(0)], where p(0) is the central density. For 
the is othermal sphere 1(0) w 2p(0)ro/Yy (Binn ev & Tremaina, 
1987) and we thus have a « 2/9 ~ 0.222. In conclusion, a a 0.2 
and the value is relatively insensitive to the choice of model 
(compared to rj). 

The central surface bright ness in the V-band is about 1(0) w 
2.5 x 10 6 L o pc" 2 dMcLaughlin & van der Marell 120051) and to- 
gether with r Q x 0.3 pc and Yy » 014 dMackev & GilmoreL 
12003b iMcLaughlin & van der MarelL 120051) we find a predicted 
central velocity dispersion of <x 1D » 5.3 km s -1 . Our measured 
dispersion is consistent with this value and we conclude that 
R136 is in virial equilibrium in the inner 5pc. This is also 
consis tent with a normal stel lar initial mass function (IMF) for 
R136 ( Andersen et al., 2009), as the expected velocity disper- 
sion would be a factor of a few lower, for example, if the IMF 
was truncated at the low-mass end (the mass-to-light ratio would 
be lower in that case). 

The expected velocity dispersion should actually be cor- 
rected for the fact that our stars are not at the centre of the cluster, 
unless we consider the cluster as an isothermal sphere in which 
case the dispersion is the same everywhere. For the Plummer 
model and the modified Hubble profile, the velocity dispersion 
for an isotropic model decreases with radius. For example, in the 
modified Hubble profile, the velocity dispersion at 5 pc should 
be ~40% lower than at the centre for a core size of 0.3 pc, so 
in virial equilibrium we would expect to measure a dispersion 
of ~3kms~ 1 at 5pc, which is still in relatively good agreement 
with our measured dispersion. 

Other effects that we have not taken into account could also 
influence our estimate of the expected velocity dispersion in 
virial equilibrium. Mass segregation, for example, would result 
in a higher central surface brightness, a smaller core size, and 
a lower mass-to-light ratio. It is however not clear what would 
be the net effect on the velocity dispersion. Radial anisotropy 
would make the (projected) velocity dispersion profile decline 
more at larger radii with respect to an isotropic model with 



1 The radius r in the isothermal model is the radius where the pro- 
jected density falls by roughly half its central value. For the EFF model 
with y = 2 the projected density is exactly half the central value at r 
and this is why we use the same symbol. 
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the sa me density profile dClarkson et al.Ll20Tll : IWilkinson et all 
2004). On the other hand, tangential anisotropy could flatten the 
profile by reducing cr 1D in the core and increasing it in the outer 
parts, which is an interesting perspective considering the evi- 
dence for rotation (Henault-Brunet et al., submitted) and the rel- 
atively flat velocity dispersion profile that we obtained for R136. 
A similarly flat profile was also found fo r the Arches cluster al - 
beit within a much smaller radial extent dClarkson et all 1201 ll) . 
Detailed numerical modeling R136 is required for a meaning- 
ful quantitative discussion of the complicated effects outlined 
above, but this is beyond the scope of the present paper. 

Other young massive clusters have recently been found 
to have a low velocity dispersion, suggesting that they are 
virial or even subvirial. Velocity dispersions of 4.5+0.8 km s~ l 
and 5.4±0.4kms~ 1 were reported, respectively, for NGC 3603 
(iRochau et all 120101) and the Arches cluster dClarkson et all 
1201 lb using proper motion measurements. From RV measure- 
ments of five yellow hypergiants and one luminous blue vari- 
able in W esterlund 1 s h owin g little RV variations over 2 to 
3 epochs, ICottaar et alj d2012l) estimated the velocity disper- 
sion of this cluster to be 2.1 + ^?kms~ 1 . From single-epoch 
near-infrared spectroscopy, Men gel & Tacconi-Garmanl (2007) 
found 5.8+2.1 km s _1 for the same cluster from the RVs of 
four red supergiants, 8.4 km s" 1 from t en post-main-sequence 
stars (M engel & Tacconi-Garmanl 120081) . and finally 9.2kms~ 1 
from a sample of four red supergi ants, five yellow hypergiants, 
and o ne B-type emission-line star dMengel & Tacconi-Garmanl 

2009) . Note that these studies of Westerlund 1 use stars of spec- 
tral types and luminosity class es that are known to be pulsators 
or int rinsic RV variables (e.g. iRitchie et all 120091 : IClark et ail 

2010) , which along with the small number statistics (both in 
terms of number of stars and number of epochs) might explain 
the range of values obtained for the velocity dispersion. 

Given the young age of these clusters for which low velocity 
dispersions were found, including R136, this might look some- 
what surprising if we expect the clusters to be expanding follow- 
ing gas expulsion. However, if the age of a cluster is at least few 
crossing times it might have had time to re-virialize, in which 
case the low velocity dispersions measured are not so surpris- 
ing. A consequence of these results is that these young massive 
clusters are certainly not being disrupted by gas expulsion, and 
in fact appear to be stable from a very young age. From now 
on, gas expulsion will not have a large effect on the dynamical 
evolution of the cluster. A gas-free cluster in virial equilibrium 
is therefore a reasonable initial condition for dynamical simula- 
tions. The gas expulsion scenario would also predict strong ra- 
dial orbits in the outer p arts of the cluster, which would result in 
a stee p decline in cr ID dClarkson et all 1201 U IWilkinson et all 
l200l . and this is not observed. Several factors could poten- 
tially explain this apparent unimportance of gas expulsion in 
earl y cluster evolution, including a high star-formation efficiency 
(e.g. lGoodwin & Bastianl [20061) the formation of these clusters 
from the merging of dyn amically cool (subvirial) substructures 
(e.g. lAllison et alll2009l) . and/or a cluster formation process re- 
sulting in a de-coupled distribution of ga s and stars that offsets 
the disruptive effec t of ga s expulsion (e.g.lFellhauer et ail [2009 ; 
iMoeckel & Clarkel 1201 1 ; Kruiiss en et all 120121) . Another con- 
sequence of these results is that even if it were true that gas ex- 
pulsion had a significant effect ~1 Myr ago, then R136 would 
have had to be incredibly dense in the past. 

R136 is often considered as an extremely dense cluster, but 
it is interesting to compare the low velocity dispersion that we 
found to the much larger line-of-sight velo city dispersion of a 
globular cluster like 47 Tucanae (11 .6 km s _1 : lMcLaughlin et all 



2006). R136 is still very young, and it will loose mass and ex- 
pand due to stellar evolution. For an adiabatic mass loss, the 
radius will grow as 1/M, so the velocity dispersion will go 
down by at least a factor of two (<x 1D oc (M/r) 1 ^ 2 oc M), but 
in reality there is g oing to be even more significant expansion 
dGieles et all 1201 Oal) which will reduce the velocity dispersion 
even more by the time R136 is as old as 47 Tucanae. 

8. Conclusions 

In an effort to determine the dynamical state of the young mas- 
sive cluster R136, we used multi-epoch spectroscopy of stars in 
the inner regions of 30 Dor. We measured RVs with a Gaussian 
fitting procedure on selected key helium lines and performed a 
quantitative assessment of the variability. Out of 41 sources for 
which spectra were extracted from the ARGUS IFU data cubes, 
16 were identified as non-variable. All of these were classified 
as O-type stars, three of which were also revealed to have com- 
posite spectra. To this sample of 16 ARGUS sources, we added 
measurements from 22 apparently single stars observed in the 
surrounding regions with Medusa-Giraffe, also as part of VFTS. 

Using this sample of 38 non-variable massive stars within 
lOpc from the centre of R 136, we computed the velocity disper- 
sion of the cluster. For the stars within 5 pc, we place an upper 
limit of 6 km s on the line-of-sight velocity dispersion. This re- 
sult does not change significantly if we exclude the few 3 <x out- 
liers, supergiant candidates, and stars having composite spectra. 
We also noted that the measured velocity dispersion of the clus- 
ter includes a small contribution of ~0.5 km s _1 from rotation. 

From Monte Carlo simulations, we established that the con- 
tribution of measurement errors to the observed velocity disper- 
sion is almost negligible. We also estimated the contribution of 
undetected binaries, which is relatively small and dominated by 
long period systems beyond the detectability range of VFTS. 
When taking errors and undetected binaries into account, we 
estimate that the true velocity dispersion of the cluster (i.e. at- 
tributable only to cluster dynamics) is between 4 and 5kms _1 
for the stars within 5 pc from the centre. 

Under basic assumptions, the expected central velocity dis- 
persion in virial equilibrium was found to be *5.3kms _1 , in 
good agreement with our measurement, and we conclude that 
R136 is in virial equilibrium. Combined with the low velocity 
dispersions found in a few other young massive clusters, our re- 
sults suggest that gas expulsion does not significantly alter the 
dynamics of these clusters. 

We would have obtained a velocity dispersion of ~25 km s _1 
if binaries had not been identified and rejected, which supports 
the suggestion that the alleged super-virial state of young star 
clusters can be expl ained by the orbital motions of binary stars 
dGieles et alll20i0bh . 
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Appendix A: Individual epochs and summary of 
radial velocity/variability analysis 



Table A.l. Individual epochs of the ARGUS observations. 



Field 


Epoch # 


Exposures 


MJD 


Al 


1 


01[a-f] 


54761.237 




2 


02 [a-f] 


54761.287 




3 


03 [g-1] 


54767.283 




4 


03 [a-f] 


54845.177 




5 


04 [a-f] 


54876.131 




6 


05 [c+d] 


55173.310 




7 


05 [a+b] 


55178.167 


A2 


1 


01 [a+b] 


54790.290 




2 


02 [a+b] 


54791.140 




3 


03 [a+b] 


54846.095 




4 


04 [a+b] 


54889.102 




5 


05 [a+b] 


55173.232 


A3 


1 


01 [a+b] 


54791.200 




2 


02 [a+b] 


54791.247 




3 


03 [a+b] 


54846.256 




4 


04 [a+b] 


54890.094 




5 


05 [a] 


55202.149 


A4 


1 


01 [a+b] 


54791.317 




2 


02 [a+b] 


54792.176 




3 


03 [a+b] 


54847.161 




4 


04 [c+d] 


54892.102 




5 


04 [e+f] 


54894.039 




6 


04 [g+h] 


54896.040 




7 


04 [a+b] 


54898.057 




8 


05 [a+b] 


55201.187 


A5 


1 


01 [a+b] 


54803.172 




2 


02 [a+b] 


54803.216 




3 


03 [a+b] 


54851.099 




4 


04 [c+d] 


54893.045 




5 


04 [a+b] 


54907.085 




6 


05 [a+b] 


55204.144 



Notes. The modified Julian date (MJD) represents the central time of 
all the exposures of a given epoch. The nomenclature of individual ex- 
posures follows that presented in Paper I. 
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Table A.2. Results of the variability tests for the ARGUS sources. 



ID 


Field 


Line 


SB 2 Epoch(s) 


1 °~ f Imax 


P(x 2 , y) 


ARV max 
[krns" 1 ] 


TVS 


Notes 


542 


A5 


He i+Hen ,14026 




12.7 


1.0000 


181.9+11.4 


y 


TVS also shows significant 

variahilitv in Rnlmpr linps nnH 

V£ll 1LIU111L V 111 Uullllt,! 1±1H,l5 tlUtl 

NIV ,14058. 






Hen ,14200 




28.9 


1.0000 


185.9+5.5 


y 


Double peak in TVS indica- 
tive of hinarv motion W/pak 
P Cygni emission. 






Hen ,14542 




36.8 


1.0000 


148.5+4.2 


y 


Double peak in TVS indica- 
tive of binary motion. Weak 
P Cygni emission. 


545 


A3 


He i+He n /14026 
Hen ,14200 
He n ,14542 




1.0 
1.1 

4.9 


0.3484 
0.5688 
1.0000 


14.4+3.0 






570 


A5 






Too low S/N for RV analysis. 


585 


A5 


He r+He n ,14026 
He ii ,14200 
He ii ,14542 
4200+4542 




0.9 
2.7 
2.5 
3.2 


0.3268 
0.9991 
0.9867 
1.0000 


32.7+9.1 






1001 


A2 














Emission-line star, no suitable 
absorption line for RV analysis. 
Variability in He n ,14542 emis- 
sion from TVS. 


1002 


A2 


He i ,14388 




2.4 


0.9136 








1003 


A2 






B [e] star, no suitable absorption 
line for RV analysis. 




A2 


He i+He n ,14026 
He 1 ,14388 
He ii ,14200 
He 11,14542 

4200+4388+4542 




2.1 
2.3 
2.4 
1.6 

3.3 


0.9600 
0.9771 
0.9196 
0.8421 
0.9986 








1005 


A2 






Too low S/N for RV analysis. 


1006 


A4 


He i+He n ,14026 
He ii ,14200 
He ii ,14542 
4200+4542 




1.2 
1.9 

3.7 
3.5 


0.0869 
0.6553 
0.9992 
0.9997 








1007 


A2 


He i+He n ,14026 
He 1 ,14388 
Hen ,14200 
He 11,14542 

4200+4388+4542 




1.9 
0.8 
2.1 
0.8 
1.4 


0.9766 
0.1244 
0.8986 
0.1502 
0.5402 








1008 


A2 


He n ,14200 
He ii ,14542 
4200+4542 




1.0 
1.2 
1.6 


0.1854 
0.4761 
0.5707 










A2 


He i+He n ,14026 
He i ,1447 1 
Hen ,14200 
Hen ,14542 
4200+454? 




0.7 
0.7 
1.6 
2.3 
2.1 


0.1257 
0.1607 
0.4610 
0.8584 
0.7858 











Ad 


Hp i+Hp tt /1 4096 
nc iTnc ii AJ-r\jz.\j 

xj ptt ,7490.0 
Up TT J4S49 
4200+4542 




1.2 
1.2 
1.0 
1.5 


0.3767 
0.2252 
0.1728 
0.4831 








1011 


A4 






Too low S/N for RV analysis. 


1012 


A2 


He i+He n ,14026 




2.0 


0.8506 












He i ,14388 




2.7 


0.9949 












He 1,14471 




1.3 


0.3517 












He ii ,14200 




1.8 


0.7524 












He ii ,14542 




1.1 


0.2058 












4200+4388+4542 




2.9 


0.9826 








1013 


A4 






Too low S/N for RV analysis. 


1014 


A2 


He i+He n ,14026 
He i ,14388 
He ii ,14200 
He ii ,14542 

4200+4388+4542 




1.2 
1.5 

2.2 
2.4 
3.0 


0.5620 
0.4265 
0.9645 
0.9268 
0.9926 









continued on next page 
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Table A.2. continued 



ID Field Line 


SB 2 Epoch(s) 


^Tmax ARVmaX TVS NOtSS 

Ikins | 


1015 A2 He i+Hen ,14026 
He U4 143 
He U4388 
Hen/14200 
Hen/14542 
4200+4542 




1.1 0.2398 
1.1 0.1787 
1.9 0.8445 

2.7 0.9953 

3.8 1.0000 69.9+16.1 
4.5 1.0000 77.7+13.9 


1016 A4 Hei+Hen.14026 
He H/14200 

Hen /14542 


? 4?, 5? 
y 2, 3, 6? 

? 6? 


1.6 0.5778 

12.9 1.0000 166.4+11.4 Y Double peak in TVS indicative 

of binary motion. 

19.2 1.0000 148.5+7.6 Y Double peak in TVS indicative 

of binary motion. 


1017 A4 Hei+Hen.14026 
He H/14200 
Hen /14542 




1.7 0.4677 

3.7 1.0000 44.0+9.1 ... Weak P Cygni emission? 
6.2 1.0000 42.1+5.7 Y Weak P Cygni emission? 


1018 A4 Hei+Hen.14026 
He H/14200 
He H/14542 
4200+4542 




1.8 0.7663 
2.1 0.5696 
1.7 0.5315 

1.9 0.8805 


1019 Al He H/14200 
Hen /14542 


y 1-4, 6, 7 
y 1-4, 6, 7 


31.6 1.0000 287.5+7.2 Y Double peak in TVS indicative 

of binary motion. 

47.4 1.0000 301.5+4.8 Y Double peak in TVS indicative 

of binary motion. 


1020 A4 Hei+Hen.14026 
He H/14200 
He H/14542 
4200+4542 




1.7 0.7254 

1.8 0.4695 

1.3 0.2682 

1.4 0.5039 


1021 A4 Hei+Hen.14026 
He 1/14471 
He H/14200 
Hen /14542 




2.8 0.9950 

3.1 0.9997 

4.7 1.0000 17.3+3.1 

5.3 1.0000 11.0+1.9 


1022 A3, A4 Hei+Hen.14026 
He H/14200 
He H/14542 




3.7 0.9994 
3.4 0.9998 
3.6 1.0000 14.2+3.2 


1023 A3 He 1/14388 
He H/14200 
He H/14542 
4200+4388+4542 




1.8 0.8739 
1.2 0.5073 
0.5 0.0229 
1.5 0.6155 


He 1/14388 
Hen/14200 
He H/14542 
4200+4388+4542 




3.1 0.9923 
2.4 0.7910 
1.3 0.3178 
2.9 0.9937 

2.2 0.9284 


1025 Al He 1/14388 




1.3 0.4262 ... ... Emission-line star. Significant 

variability in emission lines and 
P Cygni profiles from TVS. 


1026 A3 Hei+Hen.14026 
He 1/14388 
He H/14200 
He H/14542 
4200+4388+4542 




1.5 0.4289 
0.8 0.1612 
0.8 0.0743 
2.2 0.9254 
1.9 0.9063 


1027 A5 Hei+Hen.14026 
He H/14200 
Hen/14542 
4200+4542 




1.7 0.7245 
4.4 1.0000 90.3+16.4 
5.3 1.0000 74.7+12.4 
6.7 1.0000 75.2+9.9 


1028 A3 He i+He H/14026 
He H/14200 
Hen/14542 
4200+4542 




1.8 0.7310 

1.2 0.2051 
1.5 0.6237 

1.3 0.3646 


1029 A5 He i+He H/14026 
He H/14200 

Hen /14542 

4200+4542 




13.0 1.0000 113.1+7.1 y 

18.0 1.0000 97.5+4.1 Y Double peak in TVS indicative 

of binary motion. 

30.6 1.0000 102.5+2.5 Y Double peak in TVS indicative 

of binary motion. 

35.6 1.0000 101.0+2.1 



continued on next page 
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Table A.2. continued 



ID Field Line 


SB 2 Epoch(s) 


\^\ P(x\v) ARV max TVS Notes 
[kms- 1 ] 


1030 A5 




Too low S/N for RV analysis. 


1031 A5 He i+He H/14026 
He n ,14200 
Hen ,14542 


y 1-3,5-6 

Y 1-3,5-6 

✓ 1-3,5-6 


SB3 at epochs #1 and 2? 
y SB3 at epochs #1 and 2? 
3.3 1.0000 35.6+10.6 / SB3 at epochs #1 and 2? 


1032 A5 He i+Hen ,14026 
He i ,14388 
He 11,14200 
Hen ,14542 
4200+4542 




3.8 0.9995 

I. 3 0.3758 

8.5 1.0000 45.9+5.4 

II. 8 1.0000 39.3+3.4 / 
14.5 1.0000 40.7+2.9 


1033 A5 He i+Hen ,14026 
He i ,14143 
He i ,14388 
He ii ,14200 
He ii ,14542 


y 1,5 
y i 
y i 
y 1,2 
y 1,2,4 


3.4 1.0000 27.6+5.8 


1034 A5 He i+Hen ,14026 
He i ,14143 
He 1,14388 
He i ,14471 
Hen ,14200 
He ii ,14542 
4200+4388+4542 




4.7 1.0000 13.3+2.1 
3.9 0.9999 

3.7 0.9984 

5.9 1.0000 7.5+1.6 

1.8 0.7640 
3.3 0.9996 

4.9 1.0000 7.3+1.4 


Hen ,14542 
4200+4542 




1.1 0.2005 
0.8 0.0504 
1.3 0.2300 


1036 A5 




Too low S/N for RV analysis. 


1037 A5 




Too low S/N for RV analysis. 



Notes. For SB2 systems, the values presented here refer to the primary component. 



Table A.3. RVs (in kms ') for individual epochs for all the ARGUS sources 
suitable for RV analysis. 



ID 


line 


RV, 


RV 2 


RV 3 


RV 4 


RV 5 


RV 6 


RV 7 


RV g 


542 


He ii 4542 


254.1 ±2.6 


256.9 ± 3.4 


158.4 ± 2.5 


227.6 ± 3.5 


229.7 ± 2.5 


108.4 ± 2.4 






545 


He ii 4542 


230.5 ± 1.5 


230.1 ± 1.8 


228.4 ± 2.0 


232.2 ±2.1 


242.8 ± 2.2 








585 


Hen4200+Hen 4542 


280.8 ± 6.6 




248.1 ±6.3 


273.3 ± 4.5 


273.7 ± 4.0 


249.9 ± 5.3 






1002 


He i 4388 


257.6 ± 12.0 


275.2 ± 13.6 


264.9 ± 13.3 


264.6 ± 11.8 


302.1 ± 12.3 








1004 


He 1 4388+He n 4200+He n 4542 


285.7 ±4.7 


278.5 ± 5.2 


276.2 ±4.1 


258.4 ± 4.8 


273.4 ± 3.9 








1006 


Hen4200+Hen4542 


249.5 ± 9.8 


245.4 ± 11.3 


294.4 ± 8.5 


251.4 ± 13.6 


275.0 ± 16.4 


284.0 ± 17.9 


284.3 ± 13.6 


246.5 ± 8.3 


1007 


He 1 4388+He h 4200+He n 4542 


264.9 ± 2.6 


261.0 ±2.8 


262.9 ± 3.2 


263.1 ±2.9 


268.1 ± 2.8 








1008 


He ii 4200+He ii 4542 


277.3 ± 2.0 


281.0 ±2.4 


275.0 ± 2.7 


277.7 ± 2.1 


275.8 ± 1.8 








1009 


He ii 4200+He n 4542 


277.6 ± 3.0 


275.2 ± 3.3 


266.9 ± 3.5 


274.3 ± 3.7 


274.9 ± 2.5 








1010 


He ii 4200+He ii 4542 


284.2 ± 5.0 


287.1 ±6.5 


279.2 ± 5.3 


273.7 ± 7.4 


277.3 ± 7.5 


278.7 ± 7.5 


287.0 ± 6.4 


292.2 ± 6.0 


1012 


He 1 4388+He n 4200+He n 4542 


275.0 ± 6.4 


266.1 ±6.4 


267.2 ± 6.2 


247.8 ± 5.8 


269.9 ± 6.2 








1014 


He 1 4388+He n 4200+He n 4542 


269.3 ± 1.4 


266.9 ± 1.7 


268.1 ± 1.7 


266.7 ± 1.4 


262.3 ± 1.4 








1015 


He ii 4200+He ii 4542 


253.1 ±7.9 


279.7 ± 9.5 


217.2 ± 9.6 


294.9 ± 10.1 


265.4 ± 9.0 








1016 


He ii 4542 


316.9 ±4.4 


296.3 ± 5.4 


329.0 ± 4.5 


186.1 ± 6.1 


342.6 ± 5.5 


181.3 ±6.2 


329.8 ± 4.4 


320.5 ±4.1 


1017 


He ii 4542 


203.2 ± 3.3 


213.3 ±3.8 


191.7 ±3.6 


217.7 ± 3.1 


212.7 ±3.6 


227.2 ± 3.0 


202.8 ± 2.9 


185.1 ±4.3 


1018 


He ii 4200+He ii 4542 


255.4 ± 1.6 


260.9 ± 1.6 


259.8 ± 1.5 


259.0 ± 2.0 


254.6 ± 2.0 


258.2 ± 2.3 


260.0 ± 1.7 


258.1 ± 1.4 


1019 


He ii 4542 


362.2 ± 5.6 


381.9 ±6.9 


134.8 ± 4.0 


117.8 ± 4.1 


266.3 ± 3.8 


86.5 ± 3.1 


388.0 ± 3.6 




1020 


He ii 4200+He ii 4542 


265.3 ±2.1 


266.8 ± 2.8 


270.8 ± 2.5 


266.4 ± 3.0 


265.6 ± 3.1 


263.7 ± 3.4 


271.0 ±2.7 


268.0 ± 2.2 


1021 


He ii 4542 


248.6 ± 1.2 


256.8 ± 1.3 


249.8 ± 1.2 


245.8 ± 1.4 


247.2 ± 1.4 


251.2+ 1.5 


250.0 ± 1.2 


249.6 ± 1.1 


1022 


He ii 4542 (A3 pointing) 


222.2 ± 1.6 


222.5 ± 1.6 


227.0 ± 1.7 


230.4 ± 1.8 


226.0 ± 2.3 








1022 


He ii 4542 (A4 pointing) 


219.6 ±2.0 


216.2 ± 2.6 


229.6 ±2.1 


219.9 ± 2.5 


227.6 ± 2.6 


228.2 ± 2.5 


221.5 ±2.1 


217.7 ±2.2 


1023 


He 1 4388+He n 4200+He n 4542 


270.2 ± 4.4 


260.4 ± 3.9 


267.4 ± 4.7 


266.1 ±4.5 


272.0 ± 5.6 








1024 


He 1 4388+He n 4200+He n 4542 


272.3 ±6.1 


267.6 ±6.1 


245.5 ±8.1 


269.2 ±8.1 


272.6 ± 7.3 


262.8 ± 4.7 


255.6 ± 5.5 




1026 


He 1 4388+He n 4200+He ii 4542 


264.9 ± 1.7 


268.4 ± 1.8 


265.7 ±2.1 


262.3 ±2.1 


260.8 ± 2.6 








1027 


He ii 4200+He ii 4542 


292.3 ±4.2 


293.3 ± 5.9 


265.5 ± 5.2 


218.1 ±8.0 


237.9 ± 4.2 


260.4 ± 3.9 






1028 


He ii 4200+He ii 4542 


272.1 ±0.9 


272.7 ± 1.0 


271.6± 1.1 


270.3 ± 1.2 


272.1 ± 1.4 








1029 


He ii 4200+He n 4542 


290.7 ± 1.3 


292.7 ± 1.2 


227.0 ± 1.7 


287.7 ± 2.9 


262.5 ± 1.7 


328.0 ± 1.3 






1031 


He ii 4542 


274.7 ± 8.0 


278.0 ± 7.0 


242.4 ± 8.1 


256.5 ± 3.2 


251.6 ±5.3 


275.4 ± 4.0 






1032 


He ii 4200+He ii 4542 


263.9 ± 1.7 


264.2 ± 1.8 


272.7 ± 2.7 


285.6 ± 3.3 


259.9 ± 2.2 


300.6 ± 1.9 






1033 


He ii 4542 


268.3 ± 3.8 


266.2 ± 3.8 


240.7 ± 4.4 


248.8 ± 7.9 


251.2 ±4.3 


251.2 ±3.3 






1034 


He 1 4388+He n 4200+He n 4542 


263.7 ± 1.0 


263.2 ± 1.0 


263.8 ± 1.2 


270.3 ± 1.8 


264.4 ± 1.2 


270.5 ± 1.0 






1035 


He ii 4200+He ii 4542 


268.1 ±4.1 


262.7 ± 4.6 


272.1 ± 6.5 


271.3 ±9.2 


271.9 ±4.8 


269.7 ± 5.8 







Notes. The lines used for the final RV measurements of a given star are indicated. RVs for the Medusa targets are presented in Paper VIII. 
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Appendix B: Notes on individual ARGUS sources 

We comment here on selected individual sources, paying particular at- 
tention to sources with previous identifications, composite spectra, and 
also to those which appear multiple in the WFC3 image. 

- VFTS 542: This star is identified as a definite variable with a 
large amplitude from both ARGUS and Medusa observations. Its 
He ii ,14200 and He n ,14542 lines show a weak P Cygni component 
and it is classified as 02 If*/WN5 (Paper I), so even if it was not 
variable, its absolute RV could not be trusted. RV discrepancies as 
large as ~40 km s _1 are found at some epochs between Hen /i4200 
andHen/14542. 

- VFTS 545: It is also known as Mk35 dMelnickl [l985l) and classi- 
fied as 02 If*/WN5 (Paper I). There is a discrepancy of ~20 km s _I 
between Hen ,14200 and Hen ,14542, and its absolute RV also 
cannot be trusted. Low-amplitude variability is only detected in 
Hen ,14542, which is stronger and has smaller RV uncertainties 
compared to He n /14200. 

- VFTS 570: The ARGUS spectra of this source were not analysed for 
RV variability because their S/N was too low, but from the Medusa 
spectra it was found to be a definite RV variable with a large am- 
plitude (Paper VIII). Two stars appear to contribute significantly to 
the ARGUS source when comparing with the WFC3 image. 

- VFTS 585: This source was also found to be variable from the 
Medusa spectra (Paper VIII). Significant RV variability was de- 
tected from the relatively low S/N ARGUS spectra only once 
He ii ,14200 and He n ,14542 were fitted simultaneously. 

- VFTS 1001: This sourc e corre sponds to a known Wolf-Rayet 
star, R134 dFeast et all 119601) , classified as WN6(h) (e.g. 
iMassev & Hunted 1 19981) . Although it was detected as an X- 
ray source and suggest e d as a possible colliding-wind binary 
by Portegies Zwart et al. (2002), it is not known to be a binary. 
Interestingly, our TVS analysis reveals significant variability in the 
He ii ,14542 emission, but it is unclear if this is due to a normali- 
sation problem in a spectral range dominated by several emission 
lines, where the continuum is harder to define. 

- VFTS 1003: This source was found to be a new B[e]-type star in 
Paper I. The TVS analysis performed in the present work did not 
reveal any significant variability other than in the nebular emission 
lines (due to sky subtraction). 

- VFTS 1004: The WFC3 image suggests that two sources are con- 
tributing to VFTS 1004, but it does not display a composite spec- 
trum, it is not found to be variable, and the He i ,14388, He n ,14200 
and He n ,14542 lines all have consistent absolute RVs. 

- VFTS 1007: Similarly to VFTS 1004, VFTS 1007 appears multiple 
when inspecting the WFC3 image, but it is not variable, it does 
not have a composite spectrum, and Hei ,14388, Hen ,14200 and 
He ii ,14542 all have consistent absolute RVs. 

- VFTS 1014: The presence of Niv ,1405 8 and Siiv ,14089/4 116 
emission together with weak but well developed Hei singlet lines 
at 4121, 4143 and 4388 A suggests a composite spectrum (see 
Fig. IB. It . Based on the helium line diagnostics and the absence of 
Sim ,14552, the later component is identified as a mid/late O-type 
star. From the relative strength of N iv and Si iv, the other compo- 
nent is 02-4.5 (we cannot be more precise because our the ARGUS 
spectrum does not include the N v absorption region) , in agr eement 
with the 03 V classification of Massev & Hunter ( 1998). Even 
though its spectrum appears composite, this source did not show 
significant variability. The He i ,14388, He n ,14200 and He n ,14542 
lines all have consistent absolute RVs. 

- VFTS 1015: This source is clearly multiple by comparison with the 
WFC3 image and significant RV variability is found in He n ,14542. 
For some epochs, the RV of the Hei ,14388 line is clearly different 
from that of He n ,14200 and He n ,14542 lines. 

- VFTS 101 7: This source is variable. Its He n ,14200 and He n ,14542 
lines have a weak P Cygni component. A discrepancy of up to 
-30 kms" 1 is found in the RVs of Hen ,14200 and Hen ,14542 
at different epochs. 

- VFTS 1018: The presence of weak Niv ,14058 and Siiv ,141 16 
emission in combination with weak but well developed Hei sin- 
glet lines suggests a composite spectrum (see Fig. IB.lt . Based on 



the helium line diagnostics and the absence of Si m ,14552, the later 
component is identified as a mid/late O-type star.. From the rela- 
tive strength of the Niv and Siiv emission, the other component 
is classified as 02-4.5, in relatively good agree ment with the 03 
III(f*) classification of Mas sev & Hunted l[l998). Even though its 
spectrum appears composite, this source did not show significant 
variability. The Hei ,14388, Hen ,14200 and Hen /t4542 lines all 
have consistent absolute RVs. 

- VFTS 1019: This is a known high-mass binary (R13 6-038) classi- 
fied as 03 III(f*) + 08 by Massey & Hunte^ dl998l) . then revised 
as 03 V + O 6 V by IMassev et all d2002l) . The ARGUS spectra 
show obvious variability, a large RV amplitude, and SB2 profiles at 
several epochs. 

- VFTS 1022: This source corresponds to Mk37a=R136-014 
(Melni ckL 1 19851: IMassev & Hunted 119981) . classified as 04 If-H 
bv IMassev & Hunter! dl998h . but suggested as 03.5 If*/WN7 by 
ICrowth er & Walborn ( 20111) . 13 epochs (the source is on the edge 
of the A3 and A4 ARGUS pointings) made it possible to detect low- 
amplitude RV variability in He n ,14542. However, even if it had not 
been flagged as variable, this star would not have been suitable for 
our analysis of the dynamics. A discrepancy of ~15 km s -1 is found 
between the RVs of Hen ,14200 and Hen ,14542, and the RV of 
He i ,14388 is significantly larger than that of the Hen lines. 

- VFTS 1023: We classified this star as 08 III/V. IMassev & Hunted 
( 1998) classified it as 06, but at this subtype Hei+n ,14026 should 
be as deep as Hen ,14200 while Hei ,14471 should be significantly 
weaker than He n ,14542. Also, He i ,14143 and He i ,14388 should be 
much weaker than He n ,14200 and He n ,14542 respectively, which 
is not what we see. A possible exp lanation for the discrepan cy be- 
tween our classification and that of Ma ssev & Hunted (JT998) is that 
this source is an undetected single-lined spectroscopic binary. 

- VFTS 1025: This source appears multiple and the centre of the 
ARGUS position is offset between two stars in the WFC3 image, 
with the much brighter star being R 136c. It is interesting to note that 
we find significant variability in the TVS of this source (see Fig .[3>. 
R136c was identified as a probable binary ( Schnurr et al., 2009) and 
suspe cted to be a colliding-wind massive binary jcTowtheret all 
120101) . 

- VFTS 1026: When comparing with the WFC3 image, the centre of 
the ARGUS source appears offset between tw o stars. One of these is 
MH4 1. 03 IIKf*) dMassey & Hunted, Tl 998), also classified as 08: 
V by IWalborn & Blades! d!997l) . The light is probably dominated 
by MH41 (the brighter of the two stars), although we flagged VFTS 
1026 as having a composite spectrum (see Fig. IB. lb . as suggested 
by the presence of Niv /14058 and Siiv ,14089/41 16 emission to- 
gether with weak but well developed Hei singlet lines at 4121, 4143 
and 4388 A. Based on the helium line diagnostics and the absence 
of Sim ,14552, the later component is identified as a mid/late O- 
type star. From the relative strength of the Niv and Siiv emission, 
the other compo nent is classified as Q2-4 .5, in agreement with the 
classification of Mass ev & Hunted d!998h . This source is however 
not variable, and its Hei ,14388, Hen ,14200 and Hen ,14542 lines 
have consistent absolute RVs. 

- VFTS 1031: T his corresponds to R136-025 (03 V; 

IMassev & Hunted. 1 19981), which was flagged as a suspected 
variable by IMassev etal. (2002). In our ARGUS spectra, the 
Hei+Hen ,14026, Hen ,14200 and Hen ,14542 lines seem to 
display three components at some epochs. 

- VFTS 1034: This correspon ds to Mk32, which is itself a 
blend of R136-013 (Q8 111(f), IMassev & Hunten [1991 Q7.5 II 
Walbo rn & Bladesll 19971) and R136-074 (06 V IMassev & Hunted. 
1998). The variability in this source is more obvious in the 
Hei+Hen ,14026 and Hei ,14471 lines, but also significant when 
Hei ,14388, Hen ,14200 and Hen ,14542 are fitted simultaneously. 
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V. Henault-Brunet et al.: VFTS VII. A low velocity dispersion for R136 
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Fig. B.l. Non-variable ARGUS sources displaying composite spectra. Emission lines labelled are Niv /14058 and Siiv /L14089, 
41 16. Absorption lines labelled are He i ,14026, He i ,14143, He n ,14200, He i ,14388, He i ,14471, He n ,14542 and a diffuse interstellar 
band at 4428 A. 
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